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X. introduction 


The closed looop eigenstructure of a linear system has long been 
recognized as an important consideration in the design of control systems. 
Most attention in the control literature has been given to the behavior 
of closed loop eigenvalues as parameters of the control system are varied, 
and the root locus has emerged as a tool to study this behavior. The 
behavior of closed loop eigenvectors has also received a lot of attention 
because for multi-input multi-output (MXMO) systems the placement of 
eigenvectors can be used to shape transient responses, decouple inputs and 
outputs, and re jact certain types of disturbances. 

The design procedure we have in mind is to specify an asymptotic 
eigenstructure that can be achieved with high gain state feedback. Then 
the design can be simplified to choosing a single gain so that bandwidth 
constraints are satisfied. A way to do this using the linear quadratic 
regulator (LQR) is given in [1,2]. The design procedure suggested there 
is to specify an asymptotic eigenstructure, choose quadratic weights 
based on these specifications, and then vary a single parameter multi- 
plying the control weights. In this paper we propose the following 
alternative procedure: specify an asymptotic eigenstructure, choose a 
state feedback matrix based on these specifications, and then vary a 
single parameter multiplying this matrix. The alternative procedure does 
not have the advantages of guaranteed stability, phase and gain margins 
of the LQRj but is important none-the-less. Both procedures suffer the 
serious drawback of requiring full state feedback. A precursor to both 
of these procedures is in [3,4] where state feedback is used to achieve 
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a non-asymptotic eigenstructure. 

This paper starts with a raviaw of the finita and infinita zero 
structure of tha open loop system, and treats tha infinita zeros as 
• all defined quantities. The behavior of the closed loop eigenstructure 
is reviewed, and then an algorithm is presented to compute the asymptoti- 
cally infinite eigenstructure. Next we show how a MIMO system can be 
decomposed into separate single-input single-output (SISO) systems that 
have the same asymptotically infinite eigenstructures . Then we present 
the previously mentioned algorithm for achieving a desired asymptotic 
eigenstructure with full state feedback, and then we finish with several 
examples. 

Previous analysis techniques for multivariable root loci have ap- 
peared in [5,6]. in the former an algorithm is presented to compute the 
asymptotically infinite patterns of the closed loop eigenvalues. The 
first of our algorithms is similar to this, but differs importantly in 
that we compute the asymptotic behavior of the eigenvectors and we 
introduce a subspace decomposition of H n which can be used decompose the 
system into SISO parts. In the latter reference the root locus is in- 
terpreted as living on a Ri em ann surface. We have not used this ap- 
proach because of computational reasons and because we are not yet con- 
vinced tlat an engineer trying to design a control system need concern 
himself with Riemann surfaces. 

The linear systems we consider are restricted to having the same 
number of inputs and outputs, being controllable and observable, being 
nondegenerate, and having distinct finite zeros. Further restrictions are 
placed on the allowable asymptotically infinite behavior of the root locus. 


In later work we hope to remove some or all of these restrictions and to 
study In more detail the synthesis of output feedback. 

Notation 

Matrices are denoted by capital letters, scalars and vectors by lower 
T H 

case letters. A and y are the transpose of A and the Hermitian transpose 
of y. Subspaces are denoted by script letters, with the exception of 
p n . "fin A" and "ker A" are the image and kernel of the linear map A. 

The dimension of U is dim u , subspace inclusion is c > subspace inter** 
section is D , and linear combination of subspaces is U + I/. 
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II. Pole Zero Configuration of an Opm Loop System 

We are concerned with th« following linear time invariant system 


x - Ax + Bu 

(1) 

y * Cx 

(2) 


where 

x e * 
a 

u e s 

m 

ye*. 

lhe number of inputs and outputs are equal. He assume that the realiza- 
tion is minimal, which is equivalent to assuming that (A, B) is controllable 
and (A,C) is observable. The open loop system has n poles and n zeros 
associated with it, and in this section we review their dynamic inter- 
pretations. 

The open loop poles are the n eigenvalues of A. The open loop poles 
axe the complex frequencies that can appear in the output without appear- 
ing in the input. 

There are n zeros associated with the open loop system, of which 
p are finite and n-p are infinite. The zeros are sometimes called 
"transmission zeros." The finite zeros are defined to be those finite 
values of s which reduce the rank of 

[*A-sI ,1 


-C 


0 


“ 6 " 


we assume that tha finita zaroa are distinct and that tha system is not 

daganarata in tha sansa that not all values of s in tha complex plane 

reduce tha rank of tha matrix. Tha finita zeros are tha finite 
0 

solutions s^ of the generalized eigenvalue problem [7] , 


r o, *1 


On 

A-S i I B 


x i 

-C 0 



L. J 


L iJ 


i m If.*** P* (3) 


Under our assumptions the number of finite zeros is 0 <_ p £ n-m. As- 
sociated with each finite zero are right zero directions x^ and v^. it 
is also possible to define the generalized eigenvalue problem 

« 0 i»l,..., p, (4) 

and associated with each finite zero are also the left zero directions 

Oh oh 

y i and . The finite zeros are the complex frequencies "absorbed" 
by the system in the following sense. If at time t«0 the system is at 

0 o 8 - fc 

state x(0) « x^, then an input of the form u(t) - v^e 1 for t _> 0 will 
result in y(t) » 0 for t 0 [8] . 

The n-p zeroes at infinity are well defined and can be given the 
following interpretation [9] . If the input is of the form 

A ... 

u (t) - XV6 u, (t> ; (5) 

i«0 

where 5°(t) is an impulse, 6 1 (t) is a doublet, and so on; and if the 
initial state of time t - 0~ is 


[yJ H u° H ] 


A-sJl 


B 



then the output y(t) * 0 for t j> 0. The large matrix In (7) is called 
a Toeplitz matrix. The integer l is characteristic of the open loop 
system and will be interpreted in the next section. 
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HI • Thj Asymptotic Eigenatracturs of the Closed Loop Syitm 
We use output feedback of the fora 


u * - ^ Ky # (8) 

where K is an mxa invarti&ble matrix and k is a real number in the 
range 0 < k < ®, The closed loop system matrix is 

a c)1 - a - £ arc , (9) 


and its eigenvalues, right. 

and left eigenvectors 

are defined in the 

usual way by 



(A c i ■ V )x i - 0 

i*l , . . . , n 

(10) 

- V 1 ‘ 0 

i*l,...,n . 

(11) 


As k is varied from infinity down to zero the n closed loop eigenvalues 
trace out a root locus on the complex s plane and the right and left 
eigenvectors "spin" in 3t n . We are particularly interested in asymptotic 
behavior as k 0, which we now review. 

As k -► 0 p of the n branches of the root locus approach finite 
zeros. The right and left eigenvectors associated with the eigenvalues 
on the finite branches approach the right and left zero directions. 

The finite zeros and zero directions can be computed using the generalized 
eigenvalue problems (3,4), and for notational convience we group the 
solutions into 



S° ■ diag (s^, . . . , 



9- 


0 0 

[X. x 1 

x p 


- 0 0 , 


The remaining n-p branch** approach th« infinite zeros . Th* n-p 
right *ig*nv*ctor* associated with the eigenvalues on the infinite 
branches will asymptotically span the same subspace of iP spanned by all 
possible x<0“) of (6) , and a similar property can be derived for the 
left eigenvectors. 

The information we have given so far about the finite and in- 
finite eigenstructure can be computed knowing only the open loop system, 
and does not depend on the output feedback gain matrix K. Unfortunately 
this information is not enough for a good analysis of a control system, 
and this is especially true for an analysis of the asymptotically infinite 
eigenstructure. In the rest of this section we list properties of the 
asymptotically infinite eigenstructure, snst of which depend on K, and 
in the next section we show how these propertiies can be computed. 

The infinite branches of the root locus break into m patterns. The 
order of the i th pattern is n^, where “order” is defined to be the number 
of closed loop eigenvalues in the pattern. The following identity must 
be true: 


m 



and the highest order pattern is r. It turns out that r ■ l+l , where 
l is given in (5,6,7). There are n^ asymptotes in the 1 th pattern. 

The closed loop eigenvalues approach these asymptotes and have a magnitude 
approximately equal to 


10 



Tha aaymptotaa ara spacad (360/n^) • apart, aaka anglas with tha poaitiva 
raal axis of 

( arg (“*") + *360 \ * i m 0,1,..., n i »l , 

' \ / 

and hava a cantar of gravity ^ . Tha right and laft aiganvactors 

tsJi 

associatad with tha i pattam asymptotically span tha saaa subspacas 
of tt n as 

m o» n i”^ • 

ABV^f • • • i A BV^ 

T M T ao n i”^ T 00 

c r^, (AC) > (a c) • 

For notational convanianca wa group tha asymptotic props rtias into 


S* - diag(s",..., .*) 

i m 

00 09 00 

N - [V,,..., Vj 
i n 

go 00 00 

m ■ [n^,..., 






V 


( 01 , 

02 , 


11 , 


12 , 


21 

22 


9 


9 


Infill, 

[n,-l]2,.. . {n -ljm) .] 


t • • • # 
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Bach part ij of tha nulti-indax y daacritoaa tha A^BV** vactor. 

m m m « 

coaplax nuabari a^ and ^ and tha coaplax a-v actors and 
occur in coop lax conjugata pair a. 


Tha 

auat 
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IV. An Algorithm to Coapuf tha Asymptotically Infinite llgenatructura 
The problem of finding th« asymptotically infinite • i gene true tur • 

m m m 

can be concisely stated ui Given A,B,C, and K find S • H ( K t ^ and 
y. We first go through a aariaa of dafinitiona and than solva tha 
problem. Tha algorithm can ha programed on a computer using tha stable 
maaarical algorithms of BI8PAOC [10] (spacifically singular v&lua de- 
composition (SVD) and the eigenvalue problem) , but we make no claims that 
our algorithm is the "fastest” or "best". He apologize for the nota- 
tional nightmare that follows, but we do not yet know of any way to 
avoid it. 

, Define the matrices 

* CA 1*1 t • . . . r 

where r is the highest order pattern, and define the subspaces of 3R n 



1*^ * ker Gj 


• 

U ■ kar G, n...n ker G^ 

r i r 

■ »" 

V • kar G* 

* 

y - ker G* n. . .H kar G* . 

r 1 r 


The subspaces are nested so that 
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0 ■ c ... cu 0 " ** 

0 ■ V x Q ••• £ V 0 - *“ # 
and thalv dimensions ara 

■ dim ■ dim 1 * 0,'.., r . 

Dafina tha lndicaa 


■t * ‘i-i - l i 


1 * I# • • • f r f 


and It follows that 
r 


E »4 - » • 

i-1 1 

Wo rastrict our attantlon to systems for which m^ is tha nunbar of i 
ordar pattams and is tha nunbar of pattams greater than i ordar. 


th 


Naxt wa dafina tha natricas 
°i ' V i 


i"0 t 1 f • i • f r-1 


whara tha columns of form a basis for and tha columns of V i form 
a basis for Tha dimansions of both and ara nxl^. Wa dacoopoaa 
S°*, V* , and M* into 


S* ■ diag (S*, . . . , s") 

as as as 

N * [H^i ...» H^] 

oo at 00 

M ■ [M^j • • < i M^| 

as 

In general aach S. is a n i x*L block diagonal matrix, but wa rastrict our 

as 

attention to systems for which each is diagonal. Tha diagonal alemants 
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00 


>th 


ace the s^'s associated with i order patterns. Each of the and 
matrices have dim».iona and hav. a. column. the and n"'s 
associated with i“* order patterns. As we will see later. 


Oft 00 

U. ■ Im IT, + ... + Is N 
i i+l r 


V. - Da M* + ... + Im M°° 
i i+i r 


and in general 


i"0 , . . • , r-1 


00 


mu. n m n . ^ o 

i j 


Im M. n Im M. ? 0 
i D 


if*j 


In words the top relationships tell us that one of the possible sets of 

. 00 th 

basis vectors of is the v^'s associated with patterns of greater than i 

order, and that in general the v^'s are not orthogonal. A similar 

00 

statement can be made about the n ^ 1 s . The last definitions are 

t i ■ <1 K " iD i-i>' i(v r-i s i o i-i> ** r 


and the Jordan form decompositions 


[ w n "J 


w 


w 


,H 

i3 

H 


i4J 


i-l, 


The matrices « u and w. 3 have as many columns as there are 1 th order 

patterns, and W^ 2 and W i4 have as many columns as there are greater 
th 

than i order patterns . The dimensions of T . are l. , x l. . , 

* i i-l i-l 

W il and are x ra^, and W^ 2 and W^ 4 are x Finally 
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we conclude this long paragraph with the fact that If thara are no ith 
order pattema then nij-O, » U^. » 0, 

and the s”, N*, M*, w^, and W^ 3 matrices vanieh. 

The algorithm to find s** N* , H°°, y# *nd $ let 


1) use SVD to compute 0. , V. , and m. for i-1,... ,r. 

Jk a A 

2) Use » i for i*l # ...,r to compute y. 

3) Compute for i-1,..., r. 

4) Compute the Jordan form decompositions of T i for i«l, . . . ,r, 

00 09 09 

5) Compute S^i i and for r by 


* > 
S i - A i 


\ “ °i-l W il 


M i * Vl W i3 


6) Compute for m by 


'i n 


where 8, 


» H w 
Bv i 

» H n r i . 

n i a Bv i 


The following explanation of the first step may be helpful. The 


columns of form a basis for the kernel of 
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‘CB ' 

.caS. 


i*l f t • i f r*l , 


so to compute U i we can use SVD to compute an orthonoxmal basis of this 
tall skinny matrix. Likewise can be found by computing the kernel 
of 


T 

(CB) T 

* 

(CA“B) 


i*l , m m * , r-1 « 


It is not necessary to fill and with orthonormal vectors, but this 
is a convenient by-product of SVD. There are many other ways to compute 
the kernel of a matrix, and SVD is almost certainly not the method 
to use when computing the O^'s and V^'s by hand. However, there may be 
difficulty determining the dim U and SVD is the most reliable way to 
do this. If there is difficulty determining the dim U ^ then the root 
locus may have "3trange" behavior such as asymptotically infinite patterns 
that shift c liars at large radii. 

We note that we have vised the last column of a Toeplitz matrix to 
find the order of each of the asymptotically infinite root locus patterns. 
We have restricted our attention to systems for which this can be done. 

In general the entire Toeplitz matrix must be used. 

The generic case is when Rank (CB) ■ m. This should be viewed 
as a mathematical property, see [11] for a precise definition of "generic." 
If the system to be analyzed is generic then it has n-m transmission zeros 
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and m first order infinite patterns. The s"'s, v*'s , and n**s make up 
the eigenstructure of KCB. For design purposes the generic case is too 
restrictive to be of interest. This is easiest to see for SISO systems, for 
which root loci with second and higher order infinite patterns are common- 
place. 
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V. Proof of Analysis Algorithm 

In this ssction we prove the fifth step of the algorithm to find the 
asymptotically infinite eigenstructure. The proof of the sixth step we 
defer until the next section. The first four steps do not need to be 
proved. 

The proof of the first step is by induction and vises the fact that 
all of the closed loop eigenvalues s^ and the associated v_. vectors must 
satisfy 


Hex' 1 + *(s j )]v j - 0 


j - 1,..., n 


where 

<Ms) - C(sl-A) ~^B 


00 

E-r 

i-l s 



( 12 ) 


To show that (12) is true we note from (1,2,8,9,10, and 11) that 


- * 


- - 

A-S^I B 

-1 



-C -kK 


m m 

m *" 



and therefore for j»l,...,n we have that 


(A-s.I)x, +Bv - 0 
3 3 3 

CX - -kK'V - C(S.I-A)” 1 BV. 
3 

[kK _1 + Cfs.I-Aj^BlV. - 
3 3 


0 . 
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The first step of the induction proof is to show that the fifth 
step is valid for S*, N*, and M*. He assume without loss of generality 
that first order patterns exist. Equation (12) can be rewritten 


[kK 1 + ~ G. + o{r~-) ] v. ■ 0 j-1, .. . , n 

•j A 3 


As k •* 0 this becomes 


where 


[S®1 - T^V* - 0 


00 

s. - -ks . 

j j 


j-1,... ,m 


The nonzero eigenvalues of T^ correspond to first order patterns because 
the closed loop eigenvalues are solutions of 
00 


The single solution traces out a first order pattern. The right and 

00 00 

left eigenvectors of T^ are the v^'s and r^'s associated with first order 

00 

patterns, and therefore the fifth step is true for i-1. The v^'s as- 
sociated with second and higher order patterns lie in the kernel of T^, 

00 

which is Uj_. Heuristically speaking, there v^'s are not "trapped" by 
the Sj 1 term. 

OO QO 

The next step in the induction proof is to assume that , 

00 00 QO CO 

and H 1 . 1 are valid and then show that N., and H t are valid. If 

00 ao th 

N. . is valid then the v . vectors associated with > i order 

OO 

patterns form a basis for Therefore for each of these v^'s there 

exists a a such that v" . . Substituting into (12) ve get 



Multiply on tha laft to gat 


V i-1 [kK ' 1 + 4(, j )Ia i-l u j " 0 


j ■ If ... f * 


which raducaa to 


kV** k”^TJ + -i. 7^ 

i-1 i-1 i i- 



j *1 !••••£, a 
1-1 


As k+O this bscomss 


00 

(s*I — T» )w , ** 0 j*l| # # • f . 

3 i j i-1 

where 


00 

•j 



The nonzero eigenvalues of correspond to 1 th order patterns because 
the closed loop eigenvalues s^ are the i solutions of 


S j 



00 00 

The o^'s can be used to compute » U^ui^ and therefore - U i-i W ii* 

00 00 

Using a similar argument M i ■ . The v^*s associated with greater 

than i th order patterns are not "trapped" by the a”* term of (12) and 
therefore lie in U. . This completes tha proof. 
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VI. A High Sain Decomposition of the Systea into n 8 ISO Syst mm 

Thar* are a asymptotically infinite patterns of the root locus , 

and as we will now show we can decompose the KINO system into a 

SISO systems, each of which has one of the infinite patterns. After 

doing this we will prove the sixth step of the analysis algorithm 

(finding the center of gravity of the infinite patterns) . 

00 00 

Use the and vectors to define the following SISO systems i 


x ■ Ax + bjU 


y ■ 


i ■ 1,..., a 


where 


u - - y 
*i 


00 

b i - Bv i 


e i ■ \ I ' C 


lc ± - kr)^ H K" 1 v“ 


The return difference equation of each of the SISO systems is 
1 + 9 i (s)A i where 

g^e) ■ c^sl - A) 

» £ -T nT 8 CA j_1 BV® i-1,..., m. 

j-v 3 


Since the closed loop poles are the zeroes of the return difference 
equation we have that 


1 + ff i C«)/lc i - 0 
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i*l f • • • f ® 9 


which can be rewritten u 


*!<*> 


♦ 1-0 


i-1 f • • • $ A 


Carry out the long division and the result is 


i • 

i — • • ®o 


where 


1*1 f * * * f n 


(13) 


00 n i_ 00 

V» Hv 

g « -- 

1 00 V 1 00 

TljCA BV i 


• n t -1 00 

go 


V v i 


To verify that the asymptotically infinite pattern of each of the SXSO 
systems is the same as one of the patterns of the larger system we need 
only note that as k*0 (13) can be approximated by 


s” 1 - ^ 

* k *1 


i**l m 


To verify the sixth step of the analysis algorithm of section IV we 
approximate (13) by 


n i V 1 1 « 

s - f^s + £ s^ - 0 


I"! f m • # f 


m 


Tha *um of tha infinita aiganvaluas In thia pattam must aqual 0^ 
and tharafora tha cantar of gravity la 



i-1,. 




For ■ 1 tha tana "cantar of gravity" could bo mora appropriataly 
raplacod by "atartlng point." 

Sana of tha SISO ayatama nay hava conplax coaf f lcianta . Whan thia 
happana a aacond ayatam axlata which haa tha conplax conjugate coafflclanta. 
To work with only raal coafflclanta a 2x2 ayaton muat ba uaad. Tha 
cloaad loop aiganvaluaa of tha ccubinad asymptotic pattam ara solution* 
of 



1 • 

+ k *i> '* 


♦£ V 


0 , 


which 1a 







0 . 


Tha cantar of gravity of tha cotnbinad aaymptotic pattam can ba found 
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VII. An Alooritha to Synth— ii an Aawptotlc Xloanatructura with Full 
Etata Paadback 

Ht oonsidax a lyitM with full atata g aa d back t 
x - Ax ♦ Ba 

whara 

x e 

u « Jp . 

Wa iiitM that (A,B) la oontxollabla. Tha problam ia to ohooaa a faad- 
back matrix P auch that aa fc+O tha eloaad loop ayatam haa achiavad a 
daacrlJbad aiganatructusa . Oonciaaly atatad tha problan iat Givan 

0 0 oo m 

8 , X , S , H , and y find P. 

Tha aolution ia 

P - s"tO I] [X 0 , U Y P] -1 (14) 


whara 


Y 09 



00 "1 * 00 
ABV^, . . . , A BV X 

oo n 2”^ * “m" 1 n, 

ABVj , . * » t A BVj (•••> A BVj 


P » Paxsutation matrix that raorrangaa tha coluana of 

Y 

O' in an arbitary way axcapt that tha laat m columna 
of U Y ara 



n 2 -! 


<39 

'a* 
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In general tha X°, S°\ vf , and $ contain coop lax numbers. To work just 

with raal numbers than replace tha coop lax conjugate columns of x°, 

n", and 0 Y with thair raal and imaginary parts, and raplaca tha 2x2 blocks 

mm m 

m 

•i 

r i 

» « 

in 8* by 

“ » a» ■ 

to, i “ *i 

OB 00 

-B. . t R. , L 

L 

Ona way to intarprat this result is that atata faadback is usad to 
placa tha finite and infinica zeros and zero directions, and as tha 
faadback gain is increased tha cloeed loop eigenvalues and eigenvectors 
approach these zeros and zero directions. 

Tha desired asymptotic eigenetructure is not coop lately arbitrary. 

Wa restrict our attention to cases whara tha finite zeros are distinct 
and ths S matrix is diagonal. Mora fundamental are tha restrictions that 

(1) th. ccoflBx x[:, v*-» occur in 

conjugate pairs; 

0 0 

(2) for each zero direction x^ there must exist a such that 
(A-sJl)x° ♦ BvJ - 0; 

(3) the number of finite zeros must be 0 < p < n-mi 

0 Y 

(4) the columns of X and U must be linearly independent; and 

(5) the multi-index Y must specify m asymptotically infinite 
patterns whose orders total to n-p. 
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TO p r ove that (14) la trua we compute tha asymptotic eigenstructuro 
using tha previous methods daaorlbad in this papar. Onca F is fixed wa 
can treat tha stata feedb a ck problem as an output feedback problem 

X » Ax 4 Bu 


y • Fx 


u.- s xy 


Tha asymptotically finita aiganstructura must satisfy 


[: : :] [:i] 


1 ■ If ... r P 


This is trua bacausa wa haws assumed that a axists such that 
0 0 0 

(A-s i l)x 1 + »v i ■ Of and bacausa F is constructad in such a way that 
FxJ * 0. next wa show (for i-1, . . . ,m) that tha 1 th asymptotically 
infinite pattern is n^ order and has s^ and associated with it. 
By tha way F is constructad FO^P - [0 

n i“V • • • 


» 8 ) and in particular 


FA 


Bv^ - v^s^ . So if wa start with 


()tl ♦ 0(s i ))v i ■ 0 


>-1. 


i*l , . . . , n 


(15) 


$(«) ■ F(sX-A) B , and if wa 1st \> i » for 1*1,..., a. 


where 

than as k+0 (15) can be rewritten 


n. 


■ 0 


1*1 f • . . f Q 


(16) 
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Iquation (16) is true if a^ - -a" A, which it the equation for an 1 th order 
pattern. Th ii ocapletee the proof of the synthaaia algorithm. 

Ha nota that onca we hava choaan 8°, 8°, 8**, 8* # and Y thorn dona 
not exist any extra freedoa to plaoa tha cantar of gravities This 
is aasiaat to aaa for tha SZfO oaaa bacauaa it la a wall known claaaleal 
root loeus raault that tha opaa loop polaa and flnita aaroa detaraine tha 
cantar of gravity of tha aingla aayaptotioally infinite pattern. 

aquation (14) can ha tiaplifiad in tha ganario caaa (lank (CS) * m ) . 
Than wa have 0 ■ W* and (14) ba e oaaa 

F - mV (a*)" 1 10 X)[X°B)“ l . (1?) 

aquation (14) can ha simplified a van sera in tha 8180 caaa whan the 

atata apace realisation la In controllable canonical fora. For 8180 

ayataaa there la only one aayaptotically infinite pattern, it sufficaa 

to sat v" ■ l, and there does not exist any freedom to place the taro 
0 

directions x^, i-1,..., p. Tha procedure is to choose the finite zeros 
for i»l,..., p where 0 < p < n-li than fora the polyncadnal 

7T • ** * + ** ,+ s o 1 

>i 

and then tha atata feedback aatrix is 

F * IBg# ®p-l* * (18) 
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VIII . Examples 

In the first example we analyze the asymptotically infinite 
eigenstructure of a linear system. Then we use full state feedback 
to achieve a specified asymptotic eigenstructure. 

Example 1 

We use the same system as example 2 of [5] . Given the following 
A,B,C, and K matrices we use the algorithm of section IV to find the 
asymptotically infinite eigenstructure. 


12 

-60 


28 -4 


32 

36 

-54 

28 

-75 


12 30 


0 

45 

-92 

-8 

-134 


33 48 


12 

71 

0 

-8 

- 

14 


-6 21 

-18 

5 

24 










-4 


62 

- 

21 -12 

-18 

-35 

24 

64 

•a 

90 


18 24 


12 

60 

-200 

9 

- 

22 


10 -2 


12 

14 

-33 

■ 









a 



a 






6 

0 


8 


-3 

2 

3 


2 

1 


3 


11 

7 

-1 


6 

4 


0 


1 

4 

1 


0 

0 


1 

C ■ 

-9 

-11 

-1 


-3 

-2 


1 


5 

10 

2 


2 

1 


5 


—6 

- 3 

1 


2 

■ 

0 


3 


9 

- 7 

-9 

ml 
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* 1 0 0 • 

K - 0 1 0 

.0 0 1 . 

The open loop eigenvalues are -3 , - 2 , -1, 0, 1, 2, 3. To find the 
0^ and V i matrices we need 



It caui easily be verified that 




1 


1 


1 
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At this point we know there will ba a first, second, and third order 
infinite pattern. The multi-index y is therefore 

Y - (01, 02, 12, 03, 13, 23) 

where 

Dj ■ 1, n 2 ■ 2, and n 3 ■ 3 . 


The T^'s and their Jordan form decompositions are 


T x » KCB 


1 

0 

0 


-1 2 

0 0 

0 0 


1 

0 

0 


1 -1 

1 1 

0 1 J 


1 0 
0 0 
0 0 


0 

0 

0 J 


1-12 
0 1-1 
L0 0 1 


T 2 * (vjc -1 ^)" 1 vJcABU x 


r i oi 


l i o J 


l 

.0 



0 

0 . 


‘i r 
-0 - 1 . 


T 3 “ ( V 2 K ” 1 u 2 ) " 1 v 2 Ca2b0 2 * " 8 ‘ 
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Therefore we have that 



The center of gravities are 
\ - -2 
V 2 - - 1/2 
» "I • 

In summary, there is one first order pattern along the negative real 

axis with a radius of k S there is one second order pattern with 

- 1/2 

angles + 90°, radius k , and center of gravity -1/2; and there is 

- 1/2 

one third order pattern with angles 0°, ± 120° , radius 2k , and 


center of gravity -1 
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Example 2 

Consider a system with 


■ " 


p *1 

0 10 1 


0 

0 

-5 -4 0.1 1 

B - 

0 

0 

0.1 0 -11 


1 

0 

1 

in 

i 

o 

o 

o 


0 

■ft 

1 


For three different cases we specify the asymptotically infinite eigen- 
structure and then compute the state feedback matrix F. Since this is a 
generic example we use equation (17) to compute F. The asymptotically 
finite eigenstructure is the same in each case. 


0 0 

S 1.2 ■ - 3±2i x l,2 ' 

The infinite eigenstructure differs for each case . The specifications 
and the resulting F matrices are shown in Table 1. The root loci are 
shown in Figure 1. 

00 00 

In the first case we have ■ 1 and s^ ■ 2. The two first order 

00 __ 

patterns stay on the negative real axis. In the second case s 1 * l ± / 3 i 

if™ 

The two first order patterns make angles of ±120° with the positive real 
axis. In the third case s. . «» - VY ± i. Again there sure two first 

i f ^ 

order patterns, this time making angles of ±30°, and the system is un- 



stable for high gain. 
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Table 1 

Matrices Peed in Example 2 
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IX. Conclusion! 

For a restricted class of linear systems we have shown how to 
analyse the asymptotic eigenstructure as control gains get very large. 
More importantly, we have shown how to decompose MIMO systems into 
SISO systems that have the same asymptotically Infinite eigenstructures , 
and we have shown how to use state feedback to achieve a desired 
asymptotic eigenstructure . 
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